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Abstract 

The complete nucleotide sequence of plasmids pMCBFI and pMCBF6 was determined and analyzed. pMCBFI and pMCBF6 
form a novel clade within the IncP-l plasmid family designated IncP-l q. The plasmids were exogenously isolated earlier 
from a marine biofilm. pMCBFI (62 689 base pairs; bp) and pMCBF6 (66 729 bp) have identical backbones, but differ in their 
mercury resistance transposons. pMCBFI carries Tr\5053 and pMCBF6 carries Tn505S. Both are flanked by 5 bp direct 
repeats, typical of replicative transposition. Both insertions are in the vicinity of a resolvase gene in the backbone, 
supporting the idea that both transposons are "res-site hunters" that preferably insert close to and use external resolvase 
functions. The similarity of the backbones indicates recent insertion of the two transposons and the ongoing dynamics of 
plasmid evolution in marine biofilms. Both plasmids also carry the insertion sequence \SPstl, albeit without flanking repeats. 
ISPs7is located in an unusual site within the control region of the plasmid. In contrast to most known IncP-l plasmids the 
pMCBF1/pMCBF6 backbone has no insert between the replication initiation gene [trfA) and the vegetative replication origin 
(oriV). One pMCBFI /pMCBF6 block of about 2.5 kilo bases (kb) has no similarity with known sequences in the databases. 
Furthermore, insertion of three genes with similarity to the multidrug efflux pump operon mexEF and a gene from the NodT 
family of the tripartite multi-drug resistance-nodulation-division (RND) system in Pseudomonas aeruginosa was found. They 
do not seem to confer antibiotic resistance to the hosts of pMCBF1/pMCBF6, but the presence of RND on promiscuous 
plasmids may have serious implications for the spread of antibiotic multi-resistance. 
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Introduction 

Plasmids are independent from the bacterial host chromosome 
and conjugative broad-host range (BHR) plasmids can transfer 
between many different species in a bacterial community. Thus, 
plasmids are a particularly fluid part of the bacterial genome and 
one that may provide the cell with new traits that allow adaptation 
to selection forces, especially to transient environmental changes 
and challenges. Plasmids are likely to evolve differendy in different 
bacterial backgrounds [1] and might also be influenced by 
environmental conditions [2]. Thus analysis of plasmids from 
different environments is important for our understanding of these 
important mobile elements (MGE). There are a number of studies 
about plasmids in marine environments [3,4,5,6,7,8,9,10,11], as 
well as reports about the role of marine plasmids in antibiotic 
resistance [12,13,14,15] and population dynamics [16,17], for 
excellent reviews on this topic see [18,19]. Still, compared with 
many other environments, plasmids from marine habitats are not 
well studied and fully sequenced plasmids with marine origins are 
underrepresented in the databases. 

Using an exogenous isolation strategy, we have earlier isolated 
plasmids from marine bulk water, biofilms and the surface 
microlayer [20]. One of these plasmids, pMCBFI isolated from 



a biofilm, was shown to have a BHR and transferred to many 
Gram-negative bacteria, including Planctomyces maris [21]. We also 
measured comparatively high transfer rates of pMCBFI from a P. 
putida donor to indigenous bacteria, directiy in seawater [4]. 
pMCBFI was therefore an interesting BHR marine biofilm 
plasmid to characterize further. We determined the complete 
nucleotide sequences of pMCBFI and the related pMCBF6, and a 
brief summary of the plasmids was presented earlier [1]. It was 
also shown that they form a novel clade within the promiscuous 
IncP-l plasmid incompatibility group [1]. 

Since pMCBFI and pMCBF6 are the only representatives of 
their IncP-l clade, and except for pMLUAl, pMLU3 and 
pMUA4 [22], so far the only sequenced IncP-l plasmids isolated 
from marine environments, it is interesting to analyze both their 
backbone genes and their accessory genes, and to compare these 
with other IncP-l plasmids. Here we present a detailed analysis of 
pMCBFI and pMCBF6. 

Materials and Methods 

Bacterial strains, plasmids and growth conditions 

P. putida UWCl [23] and P. putida KT2440 [24] with pMCBFI 
and pMCBF6 were grown overnight at 26°C in Luria-Bertani 
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media [25] with 10 g NaCl per liter and supplemented with HgCl 2 
at 17 mg per liter. Escherichia coli with pMCBFI and pMCBF6 
were grown overnight at 37°C in the same medium supplemented 
with ampicillin (50 mg/liter). 

Mercury and antibiotic resistance tests 

To test for phenyl mercury resistance, P. putida KT2440 with 
either pMCBFI or pMCBF6 were grown overnight and spread as 
lawns on LB plates. Discs that were dipped in solutions of phenyl 
mercury chloride (saturated solution in ether) were dried and 
placed on the lawns. Plates were inoculated overnight at 30°C 
after which the diameters of clearing zones around the discs were 
recorded. P. putida KT2440 without plasmids was used as controls. 
For the antibiotic resistance tests susceptibility discs (OXOID, 
UK) were placed on bacterial lawns of E. coli HB101 [26] and 
CAG [27] with and without pMCBFI and pMCBF6 and 
incubated in 37°C overnight and clearing zones recorded as 
above. 

Sequencing and sequence analysis 

Sequencing of pMCBFI and pMCBF6 was performed earlier in 
our laboratory using standard techniques [1]. Sequence align- 
ments presented here were created using the eBioX program, the 
phylogenetic network was created using the SplitsTree program 
[28], and the similarity analysis was carried out using the SimPlot 
program [29]. The content of the genomes of the IncP-1 plasmids 
vary and all genes are not present in all plasmids [1]. The network 
constructed here was based on the largest segment that we could 
identify in which all genes were present in all analyzed plasmids 
(TraC to TraM). 

Nucleotide sequence accession numbers 

The complete sequences from pMCBFI and pMCBF6 are 
deposited with GeneBank CoreNucleotide (accession # 
AY950444 and EF 1075 16, respectively). 

Results and Discussion 

pMCBFI and pMCBF6 represent the first fully sequenced IncP- 
1 plasmids from marine environments. Generally the difference 
between the IncP-1 plasmids sequenced so far lies in the various 
insertions of resistance and catabolic genes, often as part of MGEs, 
onto a more or less conserved backbone structure. pMCBFI and 
pMCBF6 are obvious examples of this since they differ only in 
their mercury resistance (Hg R ) transposons that have inserted on 
otherwise identical backbones (Fig. 1A). Furthermore, the 
similarity of pMCBFI and pMCBF6 strongly indicates that the 
insertions of Tn 5053 in pMCBFI and Tn5058 in pMCBF6 are 
recent events. 

Plasmid backbone 

Although the focus of most previous studies has been on 
accessory genes and MGEs, it has been shown that also the IncP- 1 
plasmid backbones vary and can be divided into different 
evolutionary clades. Hitherto, twelve phylogenetic clades have 
been described for the IncP-1 plasmid group, which are designated 
a, P-l, (3-2, y, 5, s, q, k, r|, 9 [1,30,31,32,33,34,35] and two 
unnamed groups (Fig. 2). The pMCBFl/pMCBF6 plasmids 
described here are so far the only described plasmids from the 
IncP-1 <;-clade. 

As shown in Fig 2, the members of the IncP-1 plasmid backbone 
family form a star-phylogeny. Such topology is typical when there 
has been extensive recombination [36,37], exponential population 
growth [38,39], or if members of each clade have evolved in 



different distinct separated populations (i.e., long time of evolution 
since a common ancestor result in long external branches). It has 
previously been demonstrated that homologous recombination has 
contributed to the evolution of the IncP-1 backbones [1]. IncP-1 
plasmids have also been found in different environments, such as 
soil, wastewater treatment plants, and in human pathogens and 
fish pathogens [35], and as described here in a marine biofilm. 
Furthermore, it was recendy demonstrated that the IncP-1 plasmid 
backbone has evolved in, and genetically adapted to, vastiy 
different host bacterial species [1], It is thus likely that the IncP-1 
backbone at least periodically has evolved in different separated 
populations, which may explain the diversity and the star-shaped 
phylogeny of this BHR plasmid family. 

We have identified 44 open reading frames (ORFs) in the 
pMCBFl/pMGBF6 backbone. pMCBFI has 24 ORFs and 
pMCBF6 28 ORFs that are related to accessory genes or functions 
such as transposons, various resistance functions, IS elements and 
multi-drug efflux systems. There are various genetic distances 
between the pMCBFl/pMCBF6 backbone genes and the 
corresponding genes in plasmids from the other clades, which is 
probably explained by the history of recombination [1]. As an 
example, most of the backbone genes such as replication, 
conjugation and stable maintenance genes share 58-95% similar- 
ity with genes from plasmids pADP-1 [40], pB4 [41], R751 [31] 
and pTSA [42] from the (3 clade, which is the most studied IncP- 1 
clade. 

The origin of replication (oriV) of pMCBFl/pMCBF6 resembles 
other IncP- 1 plasmids. The nine TrfA DNA binding sites (iterons) 
in the vicinity of the pMCBFl/pMCBF6 oriV have the 17 bp 
consensus sequence N(T/G/G)GCCCCTCA(A/T)(A/G/C)T(A/ 
G)T(G/T)A, and have been conserved to some degree compared 
to both RK2 and R75 1 . The R75 1 iterons by comparison have the 
sequence (A/G/G)NGCCCC(T/C)C A(A/T)(G/C)(T/G)GTCA. 
The pMCBFl/pMCBF6 iterons are arranged as eight direct 
repeats and one inverted repeat, similar to the iterons close to oriV 
inR751. 

The TrfA protein is necessary for the replication of iteron 
containing theta replication plasmids (e.g. [43]) and trfA has been 
used as a base for plasmid classification [44] . We suggested earlier 
that pMCBFI and pMCBF6 had replication functions that were 
different from known plasmid incompatibility groups based on 
lack of hybridization with the inc/rep Couturier probes [20]. 
However, the general lack of similarity of the Couturier probes 
with many environmentally isolated plasmids has been noted in 
several investigations [9,18,45,46,47]. We now know that the 
IncP-1 Couturier probe only detects the IncP-1 a clade [48]. In 
fact, the nucleotide similarity between pMCBFl/pMCBF6 and 
the Couturier IncP-1 probe is only about 60% and these plasmids 
could therefore not be targeted by the probe. 

The pMCBFl/pMCBF6 backbone has no insert between the 
replication initiation gene trfA and the oriV. Except for the IncP-1 
P2 plasmids pBP136 and pAl, most other IncP-1 plasmids have an 
insert, usually several kilobases (kb), that separates trfA and oriV. 

Just as for the IncP-1 P plasmids, P MCBFl/pMCBF6 lack the 
parABCDE genes that play a key role in stability of IncP-1 a 
plasmids. Although the parABCDE genes may be absent from 
pMCBFl/pMCBF6, they apparendy possess the korA-incC-korB 
locus that mediates accurate plasmid segregation. Surface exclu- 
sion systems also protect the resident plasmid from competition, by 
preventing new plasmids to enter the cell. pMCBFl/pMCBF6 
carry the trbK homologue, which confers surface exclusion in 
plasmid RK2. We know from recombination analyses that 
pMCBFl/pMCBF6 have been involved in recombination events 
[1], which strongly suggest that these plasmids have been in the 
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Figure. 1 . Genetic map of the IncP-l q plasmids pMCBFI and pMCBF6. The backbones of the two plasmids are identical and differ only with 
regard to inserted transposons. Panel (A) shows pMCBFI with Jn5053. Coding regions are shown by arrows, indicating the direction of transcription. 
Unfilled arrows denote plasmid backbone genes, black arrows denote accessory genes. An arrow points to the insertion site of Jn5058 in the pMCBF6 
backbone. The section of the pMCBFI /pMCBF6 backbone that has a low similarity to other sequences in the databases is marked with a question 
mark. Panels B and C shows transposon Tn5053 on pMCBFI and Jn5058 on pMCBF6, respectively. Striped arrows denote genes in the transposition 
module (tni) and gray arrows denote mercury resistance (mer) genes. HP means hypothetical proteins. 
doi:10.1371/journal.pone.0092321.g001 



same host cell as other IncP-l plasmids [1]. Thus, it seems that the 
surface exclusion systems of pMCBFl/pMCBF6 have a "leaki- 
ness" that allow recombination to occur. 

One block of about 2.5 kb, between the inserted Hg 
transposons and the trb section (marked with a question mark in 
Fig. 1) has a low similarity (73% identity) to Methylophaga sp. JAM7 
[49]. 



Accessory genes, transposons and insertion sequences 

The DNA segment bearing the mercury resistance determinant 
on pMCBFI (Fig. IB) shows 99% nucleotide similarity to Tn5053 
[50]. There are five base pair direct repeats in the flanking plasmid 
DNA, which are typically formed after Tn5053 transposition [50]. 
Interestingly, Tn5053 in pMCBF 1 is inserted close to the resolvase 
gene (resA) in the backbone (Fig. 1A) which supports the suggestion 
that Tn5053 is a "resolvase-site-hunter" that insert in "hot-spots" 
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Figure. 2. Phylogenetic network of clades of the IncP-l plasmid 
family. The network is based on representative plasmids from each of 
the twelve previously described phylogenetic clades of the lncP-1 
plasmid family. The network is based on the genetic segment harboring 
the traC - traM genes for all analyzed plasmids. The two plasmids 
pMCBFI and pMCBF6 described in this study are highlighted in bold. 
Clades without designated names are marked with *. Previously 
described intra-clade recombinants [1] were not included in the 
analysis. The figure is updated from [1]. 
doi:10.1371/journal.pone.0092321.g002 

associated with res sites and resolvase genes [51]. Minakhina et al. 
suggested that the plasmid resolvase may participate in the 
transposition event, but the exact mechanism is not clear [51]. If 
so, this is an interesting case in which the transposon is dependent 
on an external site-specific recombination system. The two 
inverted repeats (IR) that form the ends of Tn5053 in pMCBFI 
have 2 and 4 mismatches compared with other published 
sequences [52]. There are also 4 mismatches between the two 
IR of Tn5053 in pMCBFI. However, this apparendy did not stop 
transposition of Tn5053 into pMCBFI. Tn5053 belongs to a 
superfamily of diverse elements including retroviruses, phage Mu, 
insertion sequences of the IS3 family, and transposons Tn552, 
Tn7 and Tn5058 [53]. They are found on plasmids within the 
IncP-l group (pPUOl, pBl 1 and pMCBFl/pMCBF6) as well as 
on other plasmid groups, such as the Enterobacter cloacae plasmid 
pELC_A [54]. 

The DNA segment bearing the mercury resistance determinant 
on pMGBF6 (Fig. 1C) is highly similar to Tn5058 in P.seudomonas 
sp. ED23-33 [55]. The length of both sequences are 12, 373 bp 
and the nucleotide similarity varies from 96 to 99% except for the 
first 257 bp where similarity is 84%, but the IR at the start of the 
element again has a perfect match with Tn5058. The element has 
25 bp inverted repeats with 2 mismatches. The complex Tn5058 
is part of the Tn5053 family (see above) [53,56]. Judging from 
comparisons with the closest relatives of Tn5058, such as Tn21, 
Tn501, Tn5053, Tn5041D, Tn5718 and others, it was suggested 
that Tn5058 was formed by a complicated series of recombination 
events involving these elements [56]. Tn5058 is also carried by 
other IncP-l plasmids such as pIJBl, pWEC911 and pTP6. 
Interestingly, Tn5058 has been detected in bacteria preserved for 
15.000-40.000 years in permafrost grounds [56]. When compared 
to several Tn5058 collected in modern days, there were high 
sequence similarities suggesting that these transposons are 



genetically stable over time. We compared the Tn5058 in plasmid 
pMCBF6 with the Tn5058 isolated from permafrost, using the 
sliding window protocol implemented in the SimPlot program. 
Our results support previous suggestions that the Tn5058 is indeed 
well conserved and there is a high nucleotide identity between the 
transposon Tn5058 isolated here and the transposon Tn5058 
isolated from permafrost in most regions. Least conserved was the 
gene tniA with a 200 bp region of only than 80% identity and merB 
and merR with regions of less than 95% identity between the 
Tn5058 in pMCBF6 and the Tn5058 isolated from permafrost. 
Tn5058 has two copies of merR, merB and merD. Interestingly, it 
was suggested that the duplications in some of the mer genes in 
Tn5041D were the result of an integration, via homologous 
recombination, of a mer containing circular DNA structure [55]. 
The circular cassette was speculated to originate from an ancestral 
donor in which the mer genes were flanked by IS elements. The 
duplication of mer genes in Tn5058 might also have been the result 
of a similar event [55]. Just as for Tn5053 in pMCBFI, Tn5058 is 
flanked by 5 base pair direct repeats in pMCBF6, which strongly 
indicate a transposition event [50,52]. These flanking direct 
repeats are also seen when Tn5058 is inserted in IncP-l plasmid 
pIJBl. The insertion of Tn5058 in pMCBF6 is only 63 bp from 
the site where Tn5053 inserted in pMCBFI, which is close to the 
resolvase gene (resA) in the plasmid backbone. The insertion site 
further confirms that Tn5058 is a "res-site hunter" [51]. 

The sequence data suggest that Tn5053 confers mercury 
resistance towards inorganic mercury while Tn5058 also carries 
resistance towards organo-mercury compounds. When the two 
plasmids were tested for phenyl mercury chloride (PhHg) 
resistance the diameter of the clearing zones around PhHg discs 
was only 6 mm for P. putida 2440(pMCBF6) but 13 mm for P. 
putida 2440(pMCBFl) and the plasmid-free P. putida 2440, 
confirming the PhHg resistance of pMCBF6 but not of pMCBF 1 . 

A 2982 bp sequence between kfrA and mexE in both pMCBFI 
and pMCBF6 shows 99% (nucleotide) similarity with insertion 
sequence ISPsl (Fig. 1A). The 24 bp IR regions with 4 mismatches 
characteristic of ISPsl [57] are found in P MCBF6/pMCBFl. 
ISPsl is related to the ISL3 family, but is larger than the usual 
ISL3 elements which are 1300-1550 bp. ISPsl has been found in 
several copies in Pseudomonas stutzeri, where it has inserted into, and 
inactivated, catabolic genes [57], and in Yersinia ruckeri plasmid 
pYRl [58]. Interestingly, ISPsl was also recendy found on 
plasmid pAMEC615 from the marine Alteromonas macleodii isolated 
from the English Channel [59]. In Pseudomonas stutzeri, ISPsl was 
flanked by eight bp direct repeats (DR), indicating a replicative 
transposition event [57], but no such flanking DR was found in 
any of pMCBFl/pMCBF6, pAMEC615 or pYRl. If DR were 
formed during a replicative transposition in these plasmids, they 
may have degenerated with time, perhaps because these sequences 
are not under selective pressure. Alternatively, the ISPviwas 
inserted by some other recombination mechanism than replicative 
transposition. The insertion point of ISPsl in pMCBFl/6, within 
the plasmid control region between kfrA and upf54.8, seems to be 
unique among IncP-l plasmids, indicating that insertion points 
other than the ones between tra and trb regions and between the 
replication region and oriV, are possible. 

A region containing a putative remnant of an insertion event is 
located between klcA and klcB. This region contains an ORF with 
closest similarity (74%) to a putative transposase, as well as an IstB 
homolog. IstB is associated with IS21 family insertion sequences. 
The function of IstB is unknown, but it may assist in transposition 
[60]. We find no DR or IR as signs of a recent insertion. Insertions 
of transposons in the klcAlklcB region are seen also in IncP-l 
plasmid RP1/RK2 [61]. 
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Adjacent to ISPsl is an inserted region that shows a high 
similarity to part of the recently sequenced plasmid pALIDE201 
from Alicycliphilus denitriftcans isolated from anaerobic sewer sludge 
[62]. Furthermore, three genes within this region show similarity 
(73-67%) to mexEF and a gene from the NodT family, which are 
part of tripartite multi-drug resistance-nodulation-division (RND) 
families. Such efflux systems are common on chromosomes of 
many Gram-negative bacteria [63], e.g. the mexEF-oprN in 
Pseudomonas aeruginosa. The IncP-1 P plasmid pB4 has a region 
with about 80% similarity to a RND efflux system in P. denitrificans, 
this was the first plasmid that was shown to carry a RND efflux 
system [41]. Adjacent to the RND region in pB4 a putative 
transposon terminus was found close to a Tn5-like transposase, 
and the authors speculated that the RND efflux system might have 
been transferred to pB4 from the chromosome of an unknown 
gram-negative bacterium by transpositional cointegrate formation 
[41]. Although ISPsl is found adjacent to the met-region in 
pMCBF6/pMCBFl, we do not know if this IS element was 
involved in the mex insertion. The MexEF-OprN system in P. 
aeruginosa confers resistance to aromatic hydrocarbons, fluoroquin- 
olones, chloramphenicol, triclosan and trimethoprim [63]. Our 
experimental data shows that when tested in disk diffusion tests for 
antibiotic sensitivity on agar plates, E. coli with pMCBFI or 
pMCBF6 did not confer a higher resistance to chloramphenicol, 
nalidixic acid or trimethoprim than plasmid free cells. Further 
studies are needed to define the possible function of MexEF-OprN 
in pMCBFl/pMCBF6. Our finding that plasmids from a marine 
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biofilm carry a multi-drug efflux system shows that pB4 is not 
unique, and that many different efflux systems may potentially be 
borne by promiscuous plasmids. The possibly serious implications 
of multi-drug efflux systems on such plasmids and the effects this 
may have on the spread of multi-resistance to human pathogens 
have been discussed [41]. 

The majority of the well-known IncP-1 plasmids, and plasmids 
in general, originate from soil environments or from a clinical 
environment and other man-made systems, such as wastewater 
treatment plants [35]. pMCBFl/pMCBF6 are among the first 
IncP-1 plasmids that are sequenced from marine environments. 
Three IncP-ls plasmids (pMLUAl, pMLUA3, pMLUA4) were 
recently isolated from the air-water interface from an estuary in 
Portugal [22,64]. Analysis of these, and other plasmids from 
marine bacteria, has often revealed new features [7,10,11,22]. 
Analysis of the pMCBFl/pMCBF6 plasmids also showed, among 
other things, that they form a new clade within the IncP- 1 plasmid 
group. 
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